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ABSTRACT 



We investigate the gas velocity dispersions of a sample of filaments recently detected as part of the Herschel Gould Belt Survey in the 
IC5 146, Aquila, and Polaris interstellar clouds. To measure these velocity dispersions, we use 13 CO, C 18 0, and N 2 H + line observations 
obtained with the IRAM 30m telescope. Correlating our velocity dispersion measurements with the filament column densities derived 
from Herschel data, we show that interstellar filaments can be divided into two regimes: thermally subcritical filaments, which have 
transonic velocity dispersions (c, & cr tot < 2c s ) independent of column density, and are gravitationally unbound; and thermally 
supercritical filaments, which have higher velocity dispersions scaling roughly as the square root of column density (cr tol oc S ' 5 ), 
and are self-gravitating. The higher velocity dispersions of supercritical filaments may not directly arise from supersonic interstellar 
turbulence but may be driven by gravitational contraction/accretion. Based on our observational results, we propose an evolutionary 
scenario whereby supercritical filaments undergo gravitational contraction and increase in mass per unit length through accretion 
of background material while remaining in rough virial balance. We further suggest that this accretion process allows supercritical 
filaments to keep their approximately constant inner widths (~ 0.1 pc) while contracting. 
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1. Introduction 



ISM: individual objects: IC5146, Aquila, Polaris - ISM: clouds - ISM: structure - evolution - sub- 



Interstellar filaments have recently received special attention, 
thanks to the high quality and dynamic range of Hersche[] 
observations (|Andre et al.||2010| |Men'shcMkov et al.|[2 010, 



Molin ari et al.||2010[ |Henning et al.||2010[ |Hill et al. 



ied with molecular line observations ( e.g., |Pineda et al!] 201 1[ 
I Hacar & Tafalla|201 l[piettinen|2012l|Bourke et al.|2012||Li &| 
Goldsmith 2012J. Herschel observations now demonstrate that 
these filaments are truly ubiquitous in MCs and have a direct link 



2011 



The submillimeter dust continuum images of nearby molecu- 
lar clo uds (MCs) taken with the SPIRE ( |Griffin et al.|2010| > and 
PACS (Poglitsch et al. 2010)cameras on board Herschel provide 
key information on both dense cores on small scales (< 0.1 pc) 
and the filamentary structure of the parent clouds on large scales 
(> 1 pc), making it possible to investigate the physical con- 
nection between these two components of the cold interstellar 
medium (ISM) and to draw a more complete picture of star for- 
mation. 

Already before Herschel, filamentary structures were known 
to be present in MCs (e.g. Schneider & Elmegre enj 1979| |Myers| 
2009) and appeared to be easily produced by any numerical 
simulation of MC evolution that includes hydrodynamic (HD) 
or magneto-hydrodynamic (MHD) turbulence (e.g., Federrath 



with the formation process of prestellar cores ( Andre et al. 2010 
|Palmeirim et al.|2013> . 

Filaments are detected with Herschel in both non- 
star-forming regions, such as the Polaris translucent cloud 



(IMen'shchikov et al.|20"T0l |Miville- Deschen es et al.|2010[|Ward- 
Thompso n et al.||2010[), and in star-forming regions ( |Konyves 
|et al. pOTO Bontemps et al. 2010| l, where they are associated 
with the presence of prestellar cores and protostars ( |Andre et al. 
|2010[ ). These findings support the view that filament formation 
precedes any star forming activity in MCs. 



|et al.|2010[[Hennebelle et al.|2008|pac Low & Klessen|2004| i. 
Individual interstellar filaments have also been recently stud 



* Based on observations carried out with the IRAM 30m telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany), and 
IGN (Spain). 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA (IPilbratt et al.|2010b. 



In a previous work (Arzoumanian et al. 201 1 1, we char- 
acterized the physical properties of the filaments detected in 
three regions observed as part of the Herschel Gould Belt sur- 
vey (HGBS - |Andre et~aT1|2010| l and found that all filaments 
share a roughly constant inner width of ~0.1 pc regardless of 
their central column density and environment. This result has 
been confirmed in other fields of the HGBS ( jPeretto et al.|2 012; 
Palm eirim et al.|2013] >. Observationally, only thermally super- 
critical filaments for which the mass per unit length exceeds the 
critical value (Mij ne > Mij necr i t )show evidence of prestellar cores 
and the presence of star formation activity, whereas thermally 
subcritical filaments (Mi; ne < Mi; ne cr jt) are generally devoid of 
Herschel prestellar cores and protostars (Andre et al. 2010). The 
critical mass per unit length, Mij ne , C rit = 2c;; /G (where c s is the 
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istothermal sound speed and G is the gravitational constant), is 
the critical value required for a filament to be gravitationally un- 
stable to radial contraction and fragmentation along its length 
(Inutsuka & Miyam afl997| l. Remarkably, this critical line mass 
Mine.crit only depends on gas temperature (Ostriker 1964). It is 
approximately equal to ~ 16 M Q /pc for gas filaments at T= 10 K, 
corresponding to c s ~ 0.2 km/s. Assuming that interstellar fila- 
ments have Gaussian radial column density profiles, an estimate 
of the mass per unit length is given by Mn ne * So x Wa where 
Wlai is the typical filament width (see Appendix A of Andre et 
al. 2010 and Arzoumanian et al. 2011) and Zo = f* m nN^ is 
the central gas surface density of the filament. For a typical fil- 
ament width of ~ 0.1 pc, the theoretical value of Mii ne cr i t ~ 
16 M /pc (for T x 10 K) corresponds to a central column den- 
sity ~ 8 X 10 21 ctrT 2 (or visual extinction Ay ~ 8). To 
first order, therefore, interstellar filaments may be divided into 
thermally supercritical and thermally subcritical filaments de- 
pending on whether their central column density iV£ is higher or 

lower than ~ 8 x 10 21 ctrT 2 , respectively. Thermally supercritical 
filaments, with average column densities > 8 x 10 21 ctrT 2 , are 
expected to be globally unstable to radial gravitational collapse 
and fragmentation into prestellar cores along their lengths. The 
prestellar cores formed in this way are themselves expected to 
collapse locally into protostars. Interestingly, the critical column 
density ~ 8 x 10 21 ctrT 2 of interstellar filaments is very close to 
the star formation threshold at a background extinction Ay ~ 6- 
9 above which the bulk of prestellar cores and young stellar ob- 



Goldsmith et al. 2008; 


Heiderman et al.|2010[|Lada et al.|2010i 


Andre etal.|2010 2011 





The results summarized above on the properties of interstel- 
lar filaments were made possible thanks to the sensitivity and 
resolution of the HGBS observations, which allowed us to detect 
structures down to Ay ~ 0. 1 and to resolve individual filaments 
in nearby regions (d<0.5 kpc), down to a typical Jeans length 
of ~ 0.02 pc. While the Herschel images already are a power- 
ful tool to assess the importance of different physical processes 
involved in the formation and evolution of interstellar filaments, 
they do not provide any information on the underlying kinemat- 
ics of the gas. 

In the present work, we complement our first results obtained 
from Herschel imaging data with molecular line observations to 
gain insight into the velocity dispersion of the interstellar fil- 
aments and reinforce our understanding of filament formation 
and evolution. We present 13 CO, C I8 0, and N2H + line observa- 
tions, taken with the IRAM 30m-telescope, toward a sample of 
filaments detected with Herschel in the IC5 146, Aquila Rift, and 
Polaris flare regions. In Sect. 2 we present the observations and 
the data reduction. Sect. 3 and 4 summarize the data analysis, 
Sect. 5 presents the results and Sect. 6 discusses their implica- 
tions in the proposed scenario for interstellar filament formation 
and evolution. 



2. IRAM observations and data reduction 

Molecular line observations of a sample of filaments in IC5146, 
Aquila, and Polaris were carried on at the IRAM-30m telescope 
at Pico Veleta, Spain, during two observing runs in March and 
August 2011. 

In March 2011, we observed the C ls O(2-l) and I3 CO(2- 
1) spectral lines with FIERA (the HEterodyne Receiver Array, 
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Fig. 1. Herschel column density map of the northern and south- 
ern streamers of IC5146 with the network of filaments high- 
lighted in blue (adapted from Fig. 3b of Arzoumanian et al. 



|201 1| ). The positions of the observed spectra are plotted as red 
squares and blue triangles for N2H + and C 18 0, respectively (the 
corresponding spectra of the positions observed for filament 2 
and filament 14 are shown in Fig.|2]i. The yellow rectangle shows 
the subcritical filament 14 which was mapped in C 18 0(2-l) and 
I3 CO(2-l);cf. Fig.||] 



Schuster et al. 2004 1 and the N 2 H + (l-0) line with EMIR (the 



Eight Mixer Receiver, a multi-band mm-wave receiver, Carter 



|et al. 2012[ ) toward the filaments of IC5 146. In addition to single- 
position spectra taken in the on-off observing mode, we mapped 
filament 14 in IC5146 (yellow rectangle in Fig.fTJ in the C 18 0(2- 
1) and 13 CO(2-l) transitions using the on-the-fly mapping mode 
with HERA. The Aquila rift and the Polaris clouds were ob- 
served in August 2011 in the C 18 O(l-0) and N 2 H + (l-0) transi- 
tions with EMIR. 

The spatial resolutions of the observations correspond to the 
12", 23" and 28" FWHM beamwidths of the 30 m antenna at 
219.6 GHz [C ls O(2-l)], 109.8 GHz [C 18 O(l-0)], and 93.2 GHz 
[N 2 H + (l-0)], respectively. We used the VESPA autocorrelator 
as backend with a spectral resolution of 20 KHz. This spectral 
resolution translates to velocity resolutions of 0.027, 0.055, and 
0.064 km/s at the frequencies of the C 18 0(2-l), C 18 O(l-0), and 
N 2 H + (l-0) transitions, respectively. 

We used the frequency-switching mode for the CO obser- 
vations and the position-switching mode for the N 2 H + observa- 
tions, using a reference position offseted by ~ 10'— 15' from each 
target position. The off position was selected from the Herschel 
images and the lack of N 2 H + emission was checked using obser- 
vations in the frequency-switching mode. During the observa- 
tions, calibration was achieved by measuring the emission from 
the sky, an ambient load and a cold load every ~ 15 minutes. The 
telescope pointing was checked and adjusted every ~ 2 hours 
and the pointing accuracy was found to be better than 3". The 
typical rms of each single position (on-off mode) spectrum is 
~0.05 K (in T^ scale). All of the data were reduced with the 
GILDAS/CLASS software package (the Grenoble Image and 
Line Data Analysis System, a software provided and actively 
developed by IRAM - http://wwwlram.fr/IRAMFR/ GILDAS). 
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Fig. 2. Spectra observed in the C ls O(2-l) and the N 2 H + (l-0) 
transitions, toward filament 14 (top) and 2 (bottom), respectively. 
The best Gaussian fits to the spectra are overlaid in blue. The 
fitted parameters are given in Table 1 . 



2.1. Filament sample 

A total number of 100 positions toward 70 filaments or seg- 
ments of filament^jwere observed in the three regions (IC5146, 
Aquila, Polaris). Most of the densest filaments, with /Vh, S 
1 x 10 22 cirT 2 , were detected in N2H + (1-0), which is a low opti- 
cal depth tracer of dense gas, as well as in C ls O and 13 CO. The 
lower column density filaments were only detected in C 18 and 
13 CO, while the faintest filaments were not detectecj^] 

The final sample discussed in the remainder of this paper 
corresponds to the 46 filaments detected in either N2H + or C 18 0. 
This sample consists of 44 individual filaments detected with 
the IRAM 30m telescope in IC5146 (11 filaments), Aquila rift 
(26 filaments), and Polaris flare (7 filaments), plus 2 filaments 



corresponding to the NGC2264C elongated clump ( |Peretto et al. 
2006 1 and the DR21 ridge (Schn eider et al.|[2010 Hennemann 
etal 20121. 



The positions of the filaments observed in IC5146 are shown 
in Fig. [T[ overlaid on the Herschel column density map. The 
fila ments obs erved in Aquil a and Polaris are shown in online 
Fig. |A.5|A.8] and Fig. |AT0] respectively. A N 2 H + (1-0) spec- 
trum and a C 18 0(2-l) spectrum toward two filaments observed 
in IC5146 are shown in Fig. [2 The positions and derived param- 
eters of all the filaments used for the analysis presented in this 



2 Due to the presence of apparent intersections of filamentary struc- 
tures in the plane of the sky, there is sometimes some ambiguity in the 
definition of an "individual filament". 

3 A few filaments observed in IC5 146 (5 filaments) and Polaris (1 fil- 
ament) were only detected in 13 CO. The velocity dispersions (~0.3 km/s 
on average) of these filaments were not included in the filament sample 
due to larger uncertainties in the optical depth of this transition. 



paper are summarized in Table 1 and the corresponding spectra 
are shown in online Appendix [A] 

3. Line-of-sight velocity dispersions 

In order to combine the velocity information obtained with vari- 
ous line tracers, we estimated the total velocity dispersion of the 
mean free particle in molecular clouds. 

We estimated the non-thermal velocity dispersion of the gas 
by subtracting the thermal velocity dispersion from the linewidth 
measured for each species, assuming that the two contributions 
are independent of each other and could be added in quadrature 
(Myers][T983[l. The thermal velocity dispersion of each species 



is given by 



0~ T (/Job.s) 



m H 



(1) 



where is the Boltzmann constant, // b s is the atomic weight of 
the observed molecule, i.e. yu b s = 29 for N 2 H + and yu b s = 30 for 
C 18 0. We adopted a gas temperature of 10 K which is the typical 
temperature of starless molecular clouds and is close to the dust 
temperature measured with Herschel for the present sample of 
filaments. The non-thermal velocity dispersion is then equal to 



0"NT = 7°"obs _ °"t^<*s) ' 



(2) 



where o- obs = AV/ V8 In 2 and AV is the measured FWHM 
linewidth of the observed spectra. The linewidths were estimated 
by single Gaussian fits to the C ls O spectra (cf. Appendix |a| for 
more details) and by multiple Gaussian fits to the seven com- 
ponents of the hyperfine multiplet of the N 2 H + transition, using 
the Gaussian HyperFine Structure (HFS) fitting routine of the 
CLASS software package. This routine derives the line optical 
depth by assuming the same excitation temperature for all hyper- 
fine components, and therefore yields an estimate of the intrinsic 
linewidth. 

The total velocity dispersion of the mean free particle of 
molecular weight fi-2.33 is finally given by: 



(3) 



where <x T (^) = ~ 0.2 km/s for T = 10 K. The typical gas 

temperature in the center of the dense filaments of our sample is 
approximately r gas ~ 10 K since the dust temperature derived 
from Herschel data is 7d ust ~ 10 K and T gas is expected to be 



well coupled to Ta ust at densities > 3 x 10 4 cirT 3 (see Galli et al. 
2002 ). As to the lower-density filaments in our sample, CO ob- 
servations directly show that their gas temperature is also close 



to 10 K (e.g., [Goldsmith et al.|2008[|Heyer et al.|2 009). 



4. Velocity dispersion along filament 14 in IC5146 

The velocity dispersions of the filaments in our sample have 
been estimated from single-position observations toward each 
filament. To assess the reliability of inferring the velocity dis- 
persion of a filament from a single position spectrum, we inves- 
tigated the variations of the velocity dispersion along the crest of 
filament 14 in IC5146 which we mapped in both 13 CO(2-l) and 
C 18 0(2-l)-seeFig.|3] 

From the 13 CO(2-l)/C 18 0(2-l) ratio map, we first estimated 
the optical depth to check whether or not the C 18 0(2-l) emission 
is optically thin. 
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Fig. 3. Left: Herschel column density map of filament 14 in IC5 146 corresponding to the yellow rectangle in Fig.[T[ in units of visual 
extinction (A v ) where Nh 2 ~ A v X 10 21 cm 2 . The resolution of the map is 36.9". The black contours correspond to column 
densities of 1 and 2 xlO 21 cm 2 . These contours are reproduced in the middle and right panels. Middle: 13 CO(2-l) integrated 
intensity map (in units of K km/s) over the LSR velocity range -2.5 to 1 km/s. The map has been smoothed to a resolution of 18.5", 
which corresponds to half the resolution of the Herschel column density map. The 13 CO(2-l) integrated emission of the filament is 
concentrated within the 10 21 crrr 2 column density contour, which traces the elongated structure of the filament. Right: Same as the 
middle panel for the C 18 0(2-l) transition. 



Assuming uniform excitation temperature along the line of 
sight, 13 CO(2-l) and C 18 0(2-l) line emissions can be expressed 
as a function of optical depth as follows: 



r:>) = ^ x [l-exp(-T>))], 



(4) 



where i corresponds to either C I8 0(2-1) or 13 CO(2-l), T*/(v) 
is the detected antenna temperature for each velocity channel u, 
T' ex is the excitation temperature which we assume to be the same 
for the two isotopomers and tL is the optical depth of the (2-1) 
transition of C ls O or 13 CO. 

The optical depth values at the systemic velocity (u sys ) of 
the cloud for each transition were estimated from the ratio of the 
13 CO(2-l) and C ls O(2-l) line intensities: 



T*/ 3 (v sys ) r e 13 [l - exp(-ri 3 (iv))] 



T*/\v sys ) CT-expt-r^xv))] 



(5) 



where r e 13 = T e 18 and t' 3 = Xt™ , with X = [ 13 CO]/[C 18 0] = 
5.5 (the assumed mean value of the abundance ratio in the local 
ISM, e.g. |Wilson&Rood|1994| ). 

Maps of the mean optical depth of filament 14 were con- 
structed from the observed integrated intensity ratio map of the 
13 CO(2-l) and C 18 0(2-l) transitions, over the velocity range 



-2.5 to +1 km/s, 



R 



13/18 



f r A 13 {v)dv / f r A 18 (v)dv 

J-2.5 J-2.5 



We estimated the mean values of r^ 8 and t^] (corresponding to 
the optical depth values integrated over the velocity range -2.5 
to +1 km/s) from the /?i3/ig ratio map by solving a velocity- 
integrated form of Eq. |5]on a pixel by pixel basis. Fig. fflshows 
the mean optical depth maps of the 13 CO(2-l) and cK)(2-l) 
transitions for filament 14 (middle and right panels, respectively) 
along with the Riyis ratio map (left panel). While the bulk of 
the C I8 emission is optically thin (t^ 8 < 0.6) along the main 
axis of the filament, the 13 CO emission is significantly optically 
thick (Tj 3 > 2) in the southern half part of the filament and only 
marginally optically thin (Tj 3 ~ 0.5 - 2) in the northern part of 
the filament. 

Fig. [5] shows the positions along filament 14 where spectra 
were extracted from the maps (left-hand side of the figure) and 
the velocity dispersions measured along the filament (right). The 
velocity dispersions derived from the 13 CO(2-l) spectra were 
corrected for the broadening effect due to the finite optical depth 
of the 13 CO(2-l) transition. 

To do so, we measured the variations of the FWHM 
linewidth, AV b s , of a simple Gaussian line model, r™ odel (u), as 
a function of optical depth. The model line spectrum had the 
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-i 7 "(d)] ) where T ex is a uniform ex- 

41n2(u-i) ) 2 



corresponds 



form r™ odel (u) = r ex [l - exp- 
ortation temperature, t(v) = tq exp 

to an intrinsically Gaussian distribution of optical depths with 
velocity, u, and AVj nt is the intrinsic FWHM linewidth of the 
model. By fitting Gaussian profiles to synthetic spectra, we es- 
timated the FWHM width AV t, s of the model as a function of 
To, and derived correction factors AVi nt / AV obs ~ 0.77 - 0.92, 
for 13 CO(2-l) optical depths ranging between 0.5 and 2.5. The 
most extreme correction would be a factor of 0.74 for a I3 CO(2- 
1) optical depth o f 3, if the [ 13 CO ]/[C 18 Q] abundance ratio has 
a value of 7 (e.g., |Wu et al.|2012 i, i.e., higher than the standard 
value in the local ISM. 

Inspection of the linewdiths measured along filament 14 
shows that the variations of the internal velocity dispersion are 
small along the crest of the filament (cf. right panel of Fig. [5}, 
with an average velocity dispersion (0.25 ± 0.03) km/s. This is 
consistent with previous observations of a few low-density fil- 
aments showing that these structures have velocity dispersions 
close to the thermal velocity dispersion ~0.2 km/s for T=10 K 
(cf. |Hily-Blant|2004l |Hily-Blant & Falgarone|2009]l and which 
do not vary much along their length (|Hacar & Taf alla 2011; 
Pineda et aLl|20TT) . Recently, |Li & Goldsmith| ( |2012| istudied 
the Taurus B213 filament and found that it is characterized by a 
coherent velocity dispersion of about ~0.3 km/s. 

These results suggest that the velocity dispersion observed 
at a single position toward a filament provides a reasonably 
good estimate of the velocity dispersion of the entire filamenf] 
Nevertheless mapping observations of a broader sample of fil- 
aments in several regions would be valuable to confirm the ro- 
bustness of this assumption. 



5. Two regimes for interstellar filaments 

In this section, we investigate the variation of internal velocity 
dispersion with column density for our sample of filaments. 

Figure [6] shows the total velocity dispersion as a function 
of central column density for the 46 filaments. Values for the 
column density and velocity dispersion of the NGC2264C elon- 
gated clump (red square in the following plots) are taken from 
|Peretto et al.| ( |2006[ ). Peretto et al. estimated the column den- 
sity from 1.2 mm dust continuum mapping observations and 
the velocity dispersion of NGC2264C from N 2 H + (l-0) observa- 
tions, both obtained with the IRAM 30m telescope. The veloc- 
ity dispersion was averaged over the entire NCG2264C clump 
and reflects the velocity dispersion of the filamentary clump as 
opposed to the embedded protostellar cores. The velocity disper- 
sion of the DR21 filament (green s quare in the following plots) is 
taken from Schneid er et al.| (|2010) who derived it from N2H + (1- 
0) observations, while the column density of DR21 comes from 
Herschel observations (|Hennemann et al. 2012)> obtained as part 



Estimated mass per unit length [M Q /pc], M line - £ W (W-O.lpc) 
1 10 100 1000 



of the HOBYS key project (Motte et al. 2010). The black tri- 
angle corresponds to filament 14 in IC5146 and the associated 
error bar comes from the dispersion of the total velocity dis- 
persion along this filament as extracted from the map shown in 
Fig. [5] The vertical error bars of the other points correspond to 
the errors in estimating the linewidth from Gaussian fits to the 
observed spectra, which are on the order of ~0.01 km/s. The er- 
ror on the measured velocity dispersion is ~ 0.03 km/s for the 

4 An interstellar filament is an elongated structure characterized by 
small column density variations along its crest (less than a factor of 3; 
see Fig.[5]-left). Therefore, the velocity dispersion variations induced by 
the trend <x lot oc E ' 5 found in Sect. 5 below for supercritical filaments 
remain small (« 2) along a given filament. 
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Fig. 6. Total velocity dispersion versus observed central column 
density (uncorrected for inclination effects): blue for IC5146, 
magenta for Aquila, cyan for Polaris, red for NGC2264C, and 
green for DR21 filaments. The squares and triangles correspond 
to velocity dispersions measured from N2H + spectra and C ls O 
spectra, respectively. The vertical error bars come from the un- 
certainties in the linewidths derived from Gaussian fits to the 
observed spectra. See text for details and uncertainties in the 
column densities. The horizontal dashed line shows the value 
of the sound speed ~ 0.2 km/s for T = 10 K. The vertical 
dashed line marks the boundary between thermally subcritical 
and thermally supercritical filaments where the estimated mass 
per unit length M\ me is approximately equal to the critical value 
Mine.crit ~ 16 M /pc for T — 10 K, equivalent to a column density 
of 8 x 10 21 crrr 2 . The grey band shows a dispersion of a factor 
of 3 around this nominal value. The dotted line running from the 
bottom left to the top right corresponds to <r tot oc Nn 2 ' 5 normal- 
ized to 0.2 km/s at the subcritical/supercritical boundary. The 
blue solid line shows the best power-law fit cr tot oc A^ ' 35 ± 014 
to the supercritical filaments. 



filaments observed in Aquila and IC5146 and and 0.1 km/s for 
DR21 and NGC2264C (representing the variation of <x tot along 
the filament crests). The errors on the filament central column 
densities come from the uncertainties in deriving column den- 
sitie^j from Herschel dust continuum maps and are typically a 
factor of ~2 mainly due to uncertain assumptions concerning the 
adopted dust opacity law (see Konyves et al. 2010 Arzoumanian 
et al. 2011 1. A recent study comparing Herschel and near-IR 
extinction data in Orion A reported a weak trend between dust 
opacity and column density, at least in the regime 1 S5 Ay % 10, 
which has been interpreted as evidence of dust grain evolution 
with environment ( Roy et al.|2013[ ). But the uncertainty induced 
by this effect remains small for most of the filaments in our sam- 
ple. 

An additional systematic effect on the column density comes 
from the random inclination angles of the filaments to the plane 
of the sky, which implies an overestimation of the intrinsic col- 



5 The column densities were derived consistently for all filaments ob- 
served with Herschel in Polaris, IC5146, Aquila, and DR21. The same 
dust opacity law was used to derive the column density of NGC2264C 
from MAMBO 1.2 mm continuum observations. The column density 
estimates for NGC2264C and DR21 are more uncertain due to the pres- 
ence of prominent cores along the filament crests (cf. Peretto et al. 2006 
Motte et al.|2007| . 
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umn densities (and masses per unit length) by ~60 % on average 
for a large sample of filaments (for more details see Appendix A 
of Arzoumanian et al. 201 1) . This systematic correction factor is 
shown as a horizontal arrow pointing to the left in Fig. [6] 

The results of our velocity dispersion measurements con- 
firm that the filaments of our sample can be divided into two 
physically distinct groups based on their central column density. 
Filaments with column densities < 8 x 10 21 cirT 2 seem to have 
roughly constant velocity dispersions close to or slightly larger 
than the sound speed (c s ^ <j tot < 2 c s ), while denser filaments 
with column densities > 8 x 10 21 cm~ 2 have velocity disper- 
sions which increase as a function of column density (cf. Fig.[6]i. 
The upper x-axis of Fig.[6]shows an approximate mass-per-unit- 
length scale derived from the bottom x-axis scale by multiplying 
the central column density of each filament by a characteristic 
filament width W n i ~ 0.1 pc (i.e., Mii ne = £oxW n i - see Sect. 1). 

It can be seen from Fig.[6]that the critical mass per unit length 
Mine.cnt ~ 16 M /pc for T ~ 10 K introduced in Sect. 1 cor- 
responds to a column density boundary which divides the fil- 
aments into two regimes where thermally subcritical filaments 
(Mii ne < M]i ne cr it) have roughly constant velocity dispersions 
with a mean value of 0.26 ± 0.05 km/s, while thermally super- 



critical filaments (My n 



Mine.crit) have velocity dispersions 



which increase as a function of projected column density as a 
power law <x tot oc A^obs ' 35 ± °' 14 - The latter relation becomes 



0"tot K N H , ,C01T 



0.41 ± 0.15 



if "intrinsic' column density values are 



used (after correcting the observed values for a ~ 60% overesti- 
mation on average due to inclination effects). The division of the 
filament sample into two regimes does not correspond to a sharp 
boundary but to a narrow border zone represented by the grey 
band in Fig. [6] which results from 1) a spread of a factor of ^ 2 
in the effective width of the filaments around the nominal value 
of 0.1 pc used to convert the critical mass per unit length to a 
critical column density, 2) uncertainties in the intrinsic column 
densities of the filaments due to a random distribution of viewing 
angles, and 3) a small dispersion in the effective gas temperature 
of the filaments, implying a small spread in the theoretical value 
of the effective critical mass per unit length (see below). 

To derive more accurate filament masses per unit length, we 
constructed radial column density profiles from the Herschel 
column density maps, taking a perpendicular cut to each filament 
at the position observed with the 30m telescope. The observed 
mass per unit length, M*, of each filament was then derived 
by integrating the measured column density profile over radius. 
The mass per unit length of DR21 corresponds to the average 
value along the crest of the filament derived by |Hennemann et aT 
( |20T2l l, while the value for NGC2264C was obtained from its 
estimated mass (1650 Mq) and length (0.8 pc). These M* es- 
timates of the mass per unit length are used in the following 
discussion instead of the simpler estimates shown at the top of 
Figure [6] 

Fig. [7k shows the total velocity dispersion of each filament 
(same y-axis as in Fig. [6]l as a function of the estimated intrinsic 
mass per unit length obtained from the observed value after ap- 
plying a 60% correction factor to account for the average over- 
estimation of the intrinsic column density due to projection ef- 
fects. This panel can be divided into two parts by comparing the 
intrinsic mass per u nit length to the virial m ass per unit length 



Mine.vir = 2o"j ot /G ( |Fiege & Pudritz 2000 1, where cr m is the 
observed total velocity dispersion (instead of the thermal sound 

The dividing line is 
= 2, where 



speed used in the expression of Mii ne ci i t ). 
defined by M\\ ns y n = 2 equivalent to an 



is the virial parameter of a filamenj^] 

Fig. [TJ) shows how the virial parameter depends on the in- 
trinsic mass per unit length for our sample of filaments. This 
plot (Fig. [TJ?) can be divided into four quadrants. The effective 
critical mass per unit length - t - 2<x^[ 2 /G divides the 

sample into subcritical filaments on the left-hand side and su- 
percritical filaments on the right-hand sid^j Conceptually, we 
define cr^j as the total velocity dispersion of a filament on the 
verge of global radial collapse. Observationally, the range of 
internal velocity dispersions cr tot measured here for thermally 
subcritical and nearly critical filaments suggests that cr^ cor- 

responds to effective temperatures T e ff = jum H o"™ /k B in the 
range 10 $ T e ff ^ 20 K. This implies an effective critical mass 
per unit length Mf? ail in the range between ~ 16 and 32 M /pc, 
which is shown as a grey band in the plots of Fig. ["/J and Fig.[7}5. 
The virial parameter, aii ne ,vir, divides the sample of filaments into 
unbound (aii ne ,vii > 2) structures in the upper part of the di- 
agram and gravitationally bound structures (aii ne ,vir < 2) in 
the lower part of the diagram. The virial parameter decreases 
as aiine.vir K ^j"™ 12 f° r the subcritical filaments, while this 
dependence is much shallower, ajine,™ K M^^ 29 , for the su- 
percritical filaments of our sample. The best-fit power law ex- 
ponent for the supercritical filaments remains almost unchanged 
(-0.28 + 0.30) when the linear fitting is performed after remov- 
ing DR21 and NGC2264C from the sample. Figure[7|3 indicates 
that thermally supercritical filaments tend to be gravitationally 
bound with aii ne ,vir = 1.0 ± 0.5 on average, and that, conversely, 
all thermally subcritical filaments are unbound with aime.vii > 2. 



6. Discussion and conclusions 

The results of our 13 CO, C I8 0, and N2H + line observations to- 
ward a sample of filaments previously detected with Herschel 
show that, to an excellent approximation, interstellar filaments 
may be divided into bound and unbound structures depend- 
ing on whether their mass per unit length Mn ne is larger or 
smaller than the thermal value of the critical mass per unit length 
Mine.crit = 2c 2 /G. Indeed, thermally supercritical filaments are 
found to be self-gravitating structures in rough virial balance 
with M* ~ Mi; ne v ; r , while thermally subcritical filaments ap- 
pear to be unbound structures with < Mn ne y u -/2 and tran- 
sonic line-of-sight velocity dispersions (c s cr tot < 2c s ). This 
confirms the usefulness of the simple gravitational instability 



criterion based on Mij ne;Cr it which was adopted by |Andre et al. 
( |20 1 0} > in their analysis of the first results from the HGBS. 

The detection, in the same cloud (e.g. Aquila), of both sub- 
critical and supercritical filaments, with transonic and supersonic 
velocity dispersions respectively, suggests that the internal ve- 
locity dispersion within a filament is not directly related to the 
level of large-scale turbulence in the parent cloud. Instead, we 
propose that the internal velocity dispersion partly reflects the 
evolutionary state of a filament (see below). 



M V i, /M where M 



Bertoldi & McKee 



19921 



«lin 



(6) 



For a sphere of uniform density a v - a 
5Ra 2 /G and M is the mass of the object (e.g 

7 In this paper, we use the terms 'supercritical' and 'subcritical' to 
characterize the mass per unit length of a filament, independently of the 
magnetic field strength (but see Fiege & Pudritz 2000 for how Mij nc cr j t 
is modified in the presence of magnetic fields). This should not be con- 
fused with the terms '(magnetically) supercritical and subcritical' used 
to characterize clouds with supercritical and subcritical mass-to-flux ra- 
tios, respectively. 
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Fig. 7. (a) Total velocity dispersion versus intrinsic mass per 
unit length (corrected for inclination effects) for the sample of 
46 filaments. The symbols and colors, as well as the y-axis of 
the plot, are the same as in Fig. [6] The horizontal segment at 
the bottom of the plot shows the typical error bar on the es- 
timated masses per unit length. Specific error bars are shown 
for NGC2264C and DR21 since they are larger 5 . The horizontal 
dashed line shows the value of the thermal sound speed ~ 0.2 
km/s for r gas = 10 K. The vertical grey band marks the the- 
oretical position of the effective critical mass per unit length 
^Unecrit ^ or effective temperatures in the range 10 < r e ff < 20 K 
(see text). The dotted line running from bottom left to top right 
corresponds to M\ ml , = Mi; ne v ; r /2, which marks the bound- 
ary between unbound and bound filaments located on the left 
and right hand side of the border line, respectively. The blue 



solid line, cr tot oc M, 



0.36 ± 0.07 



shows the best power-law fit for 



the bound filaments. The best-fit exponent remains almost un- 
changed (0.31 + 0.08) when the fitting is performed after re- 
moving DR21 and NGC2264C from the sample, (b) Corrected 
virial parameter (aime.vir) versus intrinsic mass per unit length 
for the same filaments (x-axis is the same as in panel a). The 
vertical grey band is the same as in panel a and separates sub- 
critical filaments (M°^ < Mjj? crit ) on the left from supercritical 
filaments (Mj*J > M^ ecrit ) on the right. The horizontal dashed 
line (aiine ;V ir = 2) shows the boundary between gravitationally 
unbound (aime.vh > 2) and bound (aii ne ,vir < 2) filaments, where 
aiine.vir = Mn x y n - / M°^ e . The two blue solid lines show the best 
fits for the subcritical and supercritical filaments, respectively 
(see text for details). 



6.1. An evolutionary scenario for interstellar filaments? 

In the light of the results presented in this paper, we may revisit 
the issue of the formation and evolution of interstellar filaments 
and improve the scenario proposed in our earlier work based on 
HGBS observations (e.g. Arzouma nian et al.|201 

The observed omnipresence of filamentary structures in 
gravitationally unbound complexes such as the Polaris translu- 
cent cloud (Men'shchikov et al. 2010; Miville-Deschenes et al. 



2010, War d-Thompson et al.|2 010) suggests that large-scale tur- 



bulence rather than large-scale gravity play the dominant role in 
forming interstellar filaments. The finding of a characteristic fil- 
ament width ~ 0. 1 pc ( |Arzoumanian et al.|201 1) , corresponding 
to better than a factor of ~ 2 to the sonic scale below which in- 
terstellar turbulence becomes subsonic in diffuse molecular gas 



(cf. Goodman et al. 1998 ; Falgarone et al. 2009 Federrath et al. 



2010), also supports the view that filaments may form by turbu 
lent compression of interstellar gas (Padoan et al.|200T[ ). In such 
a picture, filaments coincide with stagnation gas associated with 
regions of locally converging turbulent motions, where compres- 
sion is at a maximum and relative velocity differences at a min- 
imum, and are thus expected to have relatively low (transonic) 
internal velocity dispersions (cf. Kless en et al.|2 005 ). The line- 
of-sight velocity dispersions measured here for subcritical fila- 
ments, which all are in the range c s H cr tot % 2 c s , are consistent 
with this picture. Moreover, our 13 CO(2-l) and C ls O(2-l) map- 
ping results for filament 14 in IC5146 (cf. Fig. [3] and Fig. [5]- see 
also Hacar & Tafalla 2011 for other filaments in Taurus) sug- 
gest that subcritical filaments are velocity-coherent structures in 
MCs, with small and roughly constant velocity dispersion along 
their length (cf . |Goodman et al.|1998}|Heyer et al.|2009| for the 
definition of the concept of coherence in the context of dense 
cores). 

Since subcritical filaments are unbound, they may be ex- 
pected to widen and disperse on a turbulent crossing time. Given 
the mean velocity dispersion <x tot ~ 0.3 km/s measured here for 
subcritical filaments, the typical crossing time of a ~ 0.1-pc- 
wide filament is ~ 0.3 Myr. Subcritical filaments may be ex- 
pected to disperse on this timescale unless they are confined 



by some external pressure as suggested by Fischera & Martin 



( |2012| l. These authors modeled interstellar filaments as quasi 
equilibrium isothermal structures in pressure balanc e with a tvp 
ical ambient ISM pressure P ext ~2-5xl0 4 KcnT 3 (Fischera & 



|Martin|2012| l. The predicted widths of their model filaments are 
in rough agreement with the observed filament width of ~0. 1 pc. 
Moreover, the agreement between the expected and the observed 
widths improves when polytropic models of pressure-confined 
filaments, obeying a non-isothermal equation of state P oc p< 
with y < 1 are considered (S. Inutsuka, private communication). 

The observed increase in velocity dispersion with increas- 
ing central column density for thermally supercritical filaments, 
roughly consistent with the scaling cr tot oc TVh, ' 5 (cf. Fig.|6|, pro- 
vides interesting clues to the evolution of dense, star-forming fil- 
aments. Since the filaments on the right-hand side of Fig.|6]have 
ffiine.vir ~ 1 and are self-gravitating (see Sect. 5), they should 
be unstable to radial collapse and thus have a natural tendency 
to contract with time (cf. |Inutsuka & Miyama 1992 Kawachi 



& Hanawa|1998[ ). The shape of their radial column density pro- 
files with a high degree of symmetry about the filament major 
axis and a well-defined power law regime at large radii, ap- 
proximately consistent with a radial density structure p oc r" 2 
(jArzoumanian et al.||20"TTj |Palmeirim eFaTJ |20T3) |Hill et al.| 
2012), is consistent with a theoretical model of a nearly isother- 
mal collapsing cylinder (Kawachi & Hanawa 1998; Palmeirim 
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Fig. 8. Effective Jeans diameter (ZJj^fr = 2cr^ ot /GSo) versus fila- 
ment intrinsic central column density. The vertical dashed line 
and the grey band mark the border zone between subcritical 
and supercritical filaments (cf. Fig. [6]). The dotted line running 
from the top left to the bottom right corresponds to the ther- 
mal Jeans diameter (2Rj). The blue solid lines correspond to the 
best fits for the subcritical (Df^g cc Nn 2 * ) and supercriti- 
cal (Oj"^ °c A^h 2 °' 2 ± 01 ) filament subsamples, respectively. The 
exponent of the power-law fit for the supercritical filaments be- 
comes -0.4 + 0.1 after excluding DR21 and NGC2264C from 
the sample. In contrast to the thermal Jeans length, the effective 
Jeans diameter of the supercritical filaments remains roughly 
constant with an average value of 0.14 ± 0.07 pc. 



et al.|2013 i, and therefore also suggestive of gravitational con- 
traction. If supercritical filaments are indeed contracting, then 
their central column density is expected to increase with time (cf. 



Kawachi & Hanawa|1998| l. It is thus tempting to suggest that the 
right-hand side of Fig.|6]may, at least partly, represent an evolu- 
tionary sequence for supercritical filaments, in which central col- 
umn density corresponds to an evolutionary indicator, increasing 
with time. As the internal velocity dispersion of a filament pro- 
vides a direct measure of its virial mass per unit length and su- 
percritical filaments have My me ^ Mi; ne v i r , the trend seen on the 
right-hand side of Fig. [6] suggests that the mass per unit length 
of supercritical filaments increases with time through accretion 
of background material. More direct evidence of this accretion 
process and growth in mass per unit length for supercritical fila- 
ments exists in several cases in the form of low-density striations 
or subfilaments observed perpendicular to the main filaments 
and apparently feeding them from the side. Examples include 
the B211/B213 filament in Taurus where the typical velocities 
expected for the infalling material are consistent with the exist- 
ing kinematical constraints from CO observations (Goldsmith 
et al. 2008; Palmeirim et al. 2013), the Musca filament (Cox 
et al. 2013, in preparation), and the DR21 ridge in Cygnus X 
( |Schneideretal.|2010| i. 

If such an accretion process accompanied by an increase 
in internal velocity dispersion does indeed occur in supercriti- 
cal filaments, it may explain how these filaments can maintain 
a roughly constant inner width while contracting. In collapse 



models of nearly-isothermal cylindrical filaments (cf. Inutsuka 
|& Miyama|1992| |Kawachi & Hanawa|1998| l, the flat inner por- 
tion of the radial density profile has a radius corresponding to 
the instantaneous central Jeans length Rj ~ c^/CZo, where So is 
again the central mass column density of the model filament. In 



the presence of a nonthermal component to the velocity disper- 
sion, we may thus expect the central diameter of a contracting 
filament to be roughly given by twice the effective Jeans length, 
i.e., 2 /?j, e ff = Oj eft ~ 2 of ot /G£o. This effective Jeans diame- 
ter is plotted in Fig. [8] as a function of corrected central column 
density for the sample of filaments studied in the present paper. It 
can be seen in Fig.|8]that Dj, e ff is approximately constant for the 
thermally supercritical filaments in our sample, with an average 
value of0.14± 0.07 pc. 

In this scenario, the enhanced velocity dispersions measured 
for thermally supercritical filaments would not directly arise 
from large-scale interstellar turbulence but from the accretion of 
ambient cloud gas and the gravitational amplification of initial 
velocity fluctuations through the conversion of gravitational en- 
ergy into kinetic energy during the accretion and contraction pro- 
cess (see, e.g., the numerical simulations presented by Peretto et 
al. 2007 to explain the observed characteristics of NGC2264C). 
This would also be consistent with the concept of accretion- 
driven turbulence proposed by Klessen & Hennebelle (2010} to 
explain the origin of turbulent motions in a wide range of astro- 
physical objects. 

At the same time as they contract and accrete material from 
the background cloud, supercritical filaments are expected to 



fragment into cores and form (proto)stars ( Inutsuka & Miyama 
1997 Pon et al.|201 1 1. Indeed, for filamentary clouds, fragmen 



tation into cores and subsequent local core collapse occur faster 
than global cloud collapse ( |Kawachi & Hanawa|1998||Pon et al~] 
[20TT) [Toaia et al.||2012| >. As pointed out by |Larson| p.005) , the 
filamentary geometry is thus especially favorable for fragmen- 
tation and core formation. This is also supported by observa- 
tions of supercritical filaments which are usually found to har- 
bor several prestellar cores and Class 0/Class I protostars along 
their length (cf. Andre et al.||2010 l. Core formation is particu- 
larly well documented in the case of the B211/B213 filament in 
Taurus (Schma lzl et al.|2010||Onishi et al.|2002|l or the S erpens 
South filam ent in the Aquila Rift ( |Gutermuth et al.|2008| |Maury| 
|etal.|2011| ). 

To assess the reliability of this tentative scenario for filament 
formation and evolution, comparison with dedicated numerical 
simulations would be very valuable. More extensive molecular 
line mapping of a larger sample of filaments are also very de- 
sirable, in order to set stronger observational constraints on the 
dynamics of these structures. 
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Fig. 4. Left: Integrated intensity ratio map j^ 25 r^ 13 (u)du/ J_ 2 T^ I8 (u)dv for filament 14 in IC5146. Only those pixels where the 

column density derived from Herschel is higher thanlO 21 H2.ctrT 2 have been considered (black contour in Fig. [3]). Middle: Map 
of the mean optical depth of the 13 CO(2-l) transition. The contours in green correspond to a value of 1. Right: Map of the mean 
optical depth of the C ls O(2-l) transition. The mean optical depth of the C 18 transition is t(C 18 0) < 1 over the entire filament. 




100 -100 -200 100 -100 -200 Position along the filament ["] 

X-offset (arcsec) Y— offset fare sec) 



Fig.5. Left: Integrated 13 CO(2-l) and C ls O(2-l) intensity maps of filament 14 in IC5146 (cf. Fig. [3}. Two 13 CO(2-l) contours 
are overlaid: 1.2 and 1.5 K.km/s. The blue squares and red triangles correspond to the positions along the filament where the line- 
of-sight velocity dispersion has been measured from the 13 CO(2-l) and C I8 0(2-1) maps respectively. The two black triangles on 
the C ls O(2-l) map correspond to the positions of the central pixel of the array receiver of HERA observed in a single pointing. 
The labels indicate the offsets in arcsec relative to the centre of the map at Ra(J2000) = 21:47:56 and Dec(J2000) = 48:08:55. The 
black triangle at position (-94" , - 223") is the point which represents filament 14 in the column density vs velocity dispersion 
diagram (black triangle in figure [6). Right: Total velocity dispersions, cr tot [see Eq. ^3)], measured along the same filament. The 
x-axis goes from the south to the north of the filament. The red triangle and blue squares correspond to measurements made on the 
C 18 0(2-l) and I3 CO(2-l) maps, respectively. Individual spectra were constructed by spatially averaging the signal located within 
10" of the central position to increase the signal to noise ratio of the data. The linewidths of the 13 CO(2-l) spectra were corrected 
for the broadening due to optical depth effects (see explanations in Sect. [4j. The median value of cr tot is 0.25 km/s with a narrow 
dispersion of 0.03 km/s. 
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Table 1. Summary of derived parameters of the observed spectra toward the filaments of our sample. 



Field / tracer 


Filament 




Ra 


Dec 


N K2 


M Vmc 






rms 


°" obs 


f"lot 


(TOT' 

" tot 






(J2000) 


(J2000) 


[10 21 cm" 2 ] 


M G /pc 


[km/s] 


[K(T A *)] 


[km/s] 


[km/s] 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




1 


18 


29:39.5 


-01:50:14 


213.4 


252 


7.70 


1.43 


0.67 


0.60 


0.64 


0.01 




2 


18 


32:38.5 


-02:50:24 


40.9 


176 


7.85 


0.67 


0.22 


0.30 


0.36 


0.01 




3 


18 


33:15.2 


-02:48:14 


13.0 


35 


6.42 


0.65 


0.13 


0.11 


0.23 


0.01 




4 


18 


33:06.7 


-02:47:50 


11.6 


32 


7.04 


0.26 


0.08 


0.11 


0.23 


0.01 




5 


18 


33:22.5 


-02:42:30 


12.8 


20 


6.84 


0.16 


0.05 


0.15 


0.25 


0.01 




6 


18 


29:36.6 


-01:51:55 


13.2 


39 


7.63 


0.12 


0.05 


0.11 


0.23 


0.01 




7 


18 


29:50.4 


-01:54:47 


37.9 


123 


7.67 


1.18 


0.25 


0.18 


0.27 


0.01 


Aquila 


8 


18 


30:27.6 


-01:55:32 


35.2 


139 


7.42 


0.41 


0.13 


0.20 


0.28 


0.01 


N 2 H + (l-0) 


9 


18 


30:19.1 


-01:53:09 


33.9 


90 


7.72 


0.18 


0.11 


0.41 


0.45 


0.01 




10 


18 


29:08.9 


-01:42:40 


29.9 


121 


7.05 


0.76 


0.18 


0.21 


0.29 


0.01 




11 


18 


29:19.9 


-01:52:25 


25.1 


68 


7.27 


0.32 


0.09 


0.23 


0.31 


0.01 




12 


18 


30:03.9 


-02:04:29 


20.3 


73 


7.12 


0.68 


0.19 


0.19 


0.27 


0.06 




13 


18 


29:55.8 


-02:00:18 


74.5 


182 


7.09 


1.39 


0.60 


0.42 


0.47 


0.01 




14 


18 


29:58.1 


-01:57:38 


55.8 


265 


7.41 


0.45 


0.21 


0.44 


0.48 


0.01 




15 


18 


28:50.1 


-01:34:52 


114.1 


288 


7.90 


0.79 


0.23 


0.36 


0.41 


0.01 




16 


18 


30:03.5 


-02:02:26 


10.9 


35 


7.36 


0.37 


0.07 


0.10 


0.22 


0.01 




17 


18 


29:55.4 


-01:37:02 


12.0 


40 


8.29 


0.53 


0.12 


0.14 


0.24 


0.01 




18 


18 


29:34.6 


-01:55:15 


9.2 


38 


7.34 


0.10 


0.04 


0.15 


0.25 


0.01 




2 


21 


46:28.1 


47:17:34 


25.3 


80 


4.69 


0.61 


0.17 


0.19 


0.27 


0.01 


IC5146 


6 


21 


47:17.2 


47:33:01 


17.0 


109 


3.69 


0.32 


0.11 


0.25 


0.32 


0.01 


N 2 H + (l-0) 


12 


21 


44:58.4 


47:34:15 


14.7 


41 


4.11 


0.40 


0.13 


0.19 


0.28 


0.01 




13 


21 


47:07.9 


47:39:30 


17.1 


50 


4.35 


0.25 


0.10 


0.22 


0.30 


0.01 


Polaris 


1 


02 


00:29.0 


87:41:54 


8.0 


16 


-4.62 


0.18 


0.01 


0.13 


0.24 


0.01 


N 2 H + (l-0) 




























19 


18 


29:53.5 


-03:43:39 


7.9 


10 


4.88 


1.22 


0.06 


0.23 


0.31 


0.01 




20 


18 


28:58.2 


-01:22:51 


7.8 


28 


7.17 


1.10 


0.04 


0.30 


0.36 


0.01 




21 


18 


31:06.4 


-01:18:58 


9.9 


20 


7.77 


1.10 


0.05 


0.26 


0.33 


0.01 


Aquila 


22 


18 


33:15.8 


-02:47:15 


7.2 


51 


6.47 


1.26 


0.03 


0.13 


0.24 


0.01 


C 18 O(l-0) 


23 


18 


22:51.1 


-03:26:38 


6.3 


12 


5.78 


0.36 


0.03 


0.28 


0.34 


0.01 




24 


18 


23:31.3 


-03:03:43 


7.2 


15 


6.87 


0.89 


0.03 


0.27 


0.34 


0.01 




25 


18 


29:45.6 


-02:32:21 


7.3 


22 


6.52 


0.94 


0.03 


0.15 


0.25 


0.01 




26 


18 


32:53.1 


-01:29:42 


7.0 


21 


7.28 


0.28 


0.03 


0.16 


0.26 


0.01 




3 


21 


46:42.1 


47:44:26 


2.8 


3 


1.89 


0.80 


0.08 


0.17 


0.26 


0.01 




5 


21 


45:15.3 


47:46:52 


3.0 


11 


1.42 


0.63 


0.09 


0.24 


0.31 


0.02 


IC5146 


7 


21 


44:37.1 


47:42:39 


5.2 


11 


3.71 


0.76 


0.09 


0.29 


0.35 


0.02 


C ls O(2-l) 


10 


21 


43:29.1 


47:19:29 


4.3 


16 


4.21 


0.79 


0.06 


0.13 


0.24 


0.01 




11 


21 


43:30.9 


47:47:12 


3.1 


9 


3.29 


0.42 


0.08 


0.11 


0.23 


0.01 




1 A 

14 


21 


47:48.1 


45. U4. j i 


3 o 
i.y 


i c 
t j 


U.Z4 


0.89 


0.07 


U.lz 


0.23 


0.01 




25 


21 


44:18.8 


47:32:25 


5.7 


17 


2.32 


0.47 


0.09 


0.27 


0.34 


0.03 




2 


01 


58:59.7 


87:39:33 


9.6 


16 


-4.29 


2.65 


0.04 


0.22 


0.29 


0.01 




3 


01 


40:31.4 


87:46:19 


2.4 


15 


-4.14 


0.28 


0.04 


0.22 


0.30 


0.02 


Polaris 


4 


01 


31:21.7 


87:45:56 


2.7 


25 


-4.48 


0.27 


0.03 


0.19 


0.27 


0.02 


C 18 O(l-0) 


5 


02 


00:49.1 


87:35:18 


3.4 


15 


-4.08 


0.62 


0.04 


0.17 


0.30 


0.01 




6 


01 


50:48.0 


87:43:47 


1.5 


3 


-3.46 


0.19 


0.04 


0.19 


0.28 


0.02 




7 


03 


36:47.4 


88:08:13 


2.3 


16 


-3.39 


0.12 


0.05 


0.08 


0.26 


0.01 



Notes: Columns 3 and 4: Equatorial coordinates providing the positions of the observed spectra. Columns 5 and 6: Projected column density 
and mass per unit length for each filament measured from Herschel column density maps. The listed values have not been corrected by the 
60% systematic overestimation of true column densities and masses per unit length. All the filaments detected in N 2 H + are supercritical 
while most of the filaments observed in C 18 O(l-0) and (2-1) and not detected in N 2 H + are subcritical. Columns 7, 8 and 10: Parameters of 
single and multiple Gaussian fits to the C 18 and N 2 H + observed spectra, respectively. For the N 2 H + spectra T p!:llk (Col. 8) corresponds to 
the peak temperature of the isolated component of the N 2 H + multiplet. Column 9 gives the corresponding rms. The reference positions used 
for the sposition-switched N 2 H + observations were selected from the Herschel images and correspond to 21:47:22.6 47:45:59 for IC5146, 
01:56:30.26 87:43:28 for Polaris, 18:28:23.9 -02:00:22 for Aquila, except for filaments 3, 4 and 5 for which a reference position at 18:33:27.75 
-02:51:22 was used. Columns 10 and 11: cr obs and cr tot are derived as explained in Sect.|3] 



11 



D. Arzoumanian et al.: Formation and evolution of interstellar filaments , Online Material p 1 



Appendix A: Velocity components of the observed 
filaments 



The N 2 H + (l-0) and C 18 O(l-0)/(2-l) spectra observed toward 
the sample of filaments studied in this work are shown in 
Appendix [A] The 23 filaments detected in N2H + (l-0) in Aquila, 
IC5146, and Polaris all show a single velocity component (cf. 
Fig. A. 2 Fig. A.4JA.6 and Fig. A. 9 1. This is consistent with 
the view that the filaments are velocity coherent structures (e.g., 
|Hacar &"T afalla 201 1 ). While most of the filaments detected in 
C ls O show one velocity component, some of them show multi- 
ple (two or three) C ls O velocity components (such as filaments 
20, 22, 23, 2 5, 26 i n Aqu ila and filament 4 in Polaris - see 



Fig. |A.7|A.8l andFig. 



A.ll 



Some of these multiple C 18 veloc- 
ity components may be due to line-of-sight mixing of emissions 
not necessarily tracing the selected filament, but background 
and/or foreground structures (especially in the case of the Aquila 
Rift which is a region close to the Galactic plane). 



Assuming that interstellar filaments are velocity-coherent 
structures, the most likely velocity component associated with 
each filament was considered (selecting the velocity component 
with the velocity closest to that of neighboring filaments). Its 
properties were then used in the discussion of the results pre- 
sented in this paper (see detailed explanations in the captions of 
the figures presenting the observed spectra). 

In order to investigate how the uncertainty in selecting the 
relevant velocity component affects our results, the 6 filaments 
with spectra showing more than one C 18 velocity component 
were flagged in the sample. Fig. |A.l| shows the total velocity dis- 
persion as a function of central column density for the resulting 
subsample of 38 filaments, where the velocity dispersions of the 
filaments inferred from the C ls O spectra (which show multiple 
velocity components) corresponding to filaments 20, 22, 23, 25, 
26 in Aquila and 4 in Polaris are flagged in yellow. Our general 
conclusions on filament properties, such as total velocity disper- 
sion, central column density, and mass per unit length remain 
unchanged whether or not we take the 6 flagged filaments into 
account. 



Estimated mass per unit length [M Q /pc], M line = £ W (W=0.1pc) 
10 100 



0.4 - 



0.2 - 



■ IC5146 






Aquila 






_ Polaris 




• 

i * * 


{ i 5 = 


! j: 


* t • 


sound speed 







10 10 10_ 

Central column density, N H [cm ] 

Fig.A.l. Total velocity dispersion versus observed central col- 
umn density. Similar plot to Fig.|6]for the filaments observed in 
IC5146, Aquila and Polaris. The six C 18 spectra showing two 
or three velocity components are flagged in yellow. 
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Fig. A.2. N2H + (l-0) spectra observed toward filaments 6, 12, 13 in IC5146. The filament numbers are indicated at the upper right 
of each panel. The corresponding one velocity hyperfine structure Gaussian fits are highlighted in blue. 
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Fig. A.3. C 0(2-1) spectra observed toward 6 filaments in IC5146. The filament numbers are indicated at the upper right of each 
panel. The corresponding single-component Gaussian fits are highlighted in blue. 
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Fig.A.4. N2H + (l-0) spectra observed toward 9 filaments in Aquila. The numbers indicated on the upper right hand side of the 
plots, correspond to the filaments marked on the column density map of Fig. A. 5 and listed in Table 1 . The corresponding single- 



component hyperfine structure Gaussian fits are highlighted in blue. Note that the y-axis scale is not the same for all spectra, but has 
been adjusted to match the peak temperatures which vary between ~ 0.2 to > 3 K. 
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Fig. A.5. Curvelet component of the column density map of the Aquila field around the Serpens South filaments take n from | Andre 
( 201Q| >. The curvelet component is obtained using a morphological component analysis algorithm (MCA, from Starck et al. 
i which enhances the contrast of the filamentary structure against the more diffuse background of the cloud. The positions of 



2003 



the observed spectra are plotted in red squares and blue triangles for N2H + and C I8 0, respectively. The numbers coiTespond to the 
filaments listed in Table 1 and are given at the upper right of the spectra in Fig. A.4 and some of the spectra of Fig. A. 7 
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Fig. A.6. Continuation of Fig. A.4 presenting the N2H + (l-0) spectra observed toward 9 filaments in Aquila. 
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Fig. A.7. C 18 O(l-0) spectra observed toward 6 filaments in Aquila. The corresponding Gaussian fits are highlighted in blue. Note 
that the scale of the y-axis has been adjusted to match the peak temperature of each spectrum. Note 1: Filament 20 shows two 
components separated by 0.63 km/s. The component at Vlsr = 7.17 km/s have been chosen to be tracing the filament, since this 
component has a velocity component which is closer to the average velocity (Vlsr = 7.4 km/s) of the filaments in Aquila which 
show one velocity component. Note 2: The spectrum observed at the position of filament 22 shows two velocity components, one 
centered at 6.5 km/s and the other at 7.7 km/s. The strongest component has been selected for similar reasons to that explained for 
filament 20. Note 3: The spectrum observed at the position of filament 25 has two velocity components (at 5.0 km/s and 6.52 km/s, 
respectively). The strongest component has been selected for similar reasons to that explained for filament 20. Note 4: The spectrum 
observed at the position of filament 26 shows three velocity components centered at 7.28 km/s, 8.45 km/s and 9.17 km/s, respectively. 
The component centered at 7.28 km/s has been selected because it is the component with the closest Vlsr value to the mean Vlsr 
of neighboring filaments in Aquila. 
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Fig.A.8. C 18 O(l-0) spectra observed toward filaments 23 and 24 in Aquila, located in the south west part of the region. The 
positions of the two filaments are shown on a blow up of the column density map (cf., Kon yves et al.|2010] l on the right hand side. 
Filament 24 shows a single velocity C 18 O(l-0) spectrum at Vlsr = 6.87 km/s, while C 15 O(l-0) spectrum observed toward filament 
23 shows three velocity components at Vlsr = 5.77 km/s, 9.13 km/s and 10.3 km/s. The relevant velocity component associated to 
filament 23 is most probably the component at Vlsr = 5.77 km/s (the closest to the LSR velocity of the neighboring filament 24), 
while the other components are likely associated to clouds/filaments observed on the same line-of-sight. 
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Fig.A.9. N2H + (l-0) spectrum observed toward filament 1 in Polaris. The corresponding hyperfine structure Gaussian fit is high- 
lighted in blue. 
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Fig. A. 10. Curvelet component of the column density map of a subregion in the Polaris Cloud taken from 
The positions of the observed spectra are plotted in red squares and blue triangles for N2H + and C ls O, respectively. The numbers 
correspond to the filaments listed in Table 1. 
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Fig. A.ll. C 18 O(l-0) spectra observed toward 6 filaments in Polaris. The corresponding Gaussian fits are highlighted in blue. Note 
1: The spectrum observed toward filament 4 shows two velocity components. Both velocity components have similar linewidths 
(0.19 km/s and 0.14 km/s, respectively), they are separated by ~0.8 km/s. They most probably belong to two filaments which are 
seen on the same line-of-sight (but a large map would be needed to study the kinematics of the field). The component observed at 
-4.5 km/s has been taken to be representative of filament 4. 



